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Abstract 

We present a quenched calculation of lowest bottomonium hybrid states on 
an anisotropic lattice with Landau mean-link tadpole improvement, using the im- 
proved NRQCD Hamiltonian. We investigate the quark-spin s = 0, 1 sectors which 
contain both the non-exotic 1 and the exotic 1 and demonstrate their degen- 
eracy in the case of the lowest order Hamiltonian. Both states are found at around 
1.6 GeV above the T ground state. We examine the spin-splitting for several hy- 
brid states (l -- , 1 '", 2 h ) which is due to the —c\cr.B/2mb term in the 
NRQCD Hamiltonian, all other terms having negligible effect. The spin splittings 
are well resolved outside errors and are surprisingly large. We investigate their 
dependence on c\ for several values of c\. We calculate one contribution to these 
splittings using the Born-Oppenheimer picture for hybrids and show that the ob- 
served size of the effect is plausibly explained by mixing with hybrid states with 
more than one constituent gluon. 
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1 Introduction 



The spectrum of hybrid mesons provides fresh insight into the structure of QCD. Hybrid 
mesons are therefore of intense interest, both experimentally and theoretically. In this 
paper we study heavy hybrids using lattice QCD. The results invite testing in future 
£>-factory experiments. 

Because of the large mass of the 6-quark we can apply the NRQCD approach, which 
has been successful for the T-system ]||, with some confidence. Indeed we expect 
the quarks in heavy hybrid states to be even more slowly moving than is the case for 
the T itself. We work in the quenched approximation for the gauge fields. However it 
is unlikely that the use of dynamical gauge configurations would radically change our 
results. 

In order to reconcile computational efficiency with the need for accurate measure- 
ment of the relatively large values of hybrid excitations, we use an anisotropic lattice 
that is spatially coarse but refined in the time direction. The gauge field configurations 
are determined by the spatially improved gluon field action 0| , 

r, _ _n \ -l/5 Pg iS > 1 R SS ,s' 1 Rs's',s \ [4 Pg£ 1 Rss,t 

n ~ P Uui 12 12 „6 J ^/U^uf 12u}ti? 

(1.1) 

where s, s' run over spatial links, P S;S / and P s> t are lxl plaquettes, R S s,s> and R S s,t are 
2x1 plaquettes, x is t ne bare anisotropy. 

The tadpole parameters u s and ut are defined to be respectively, the expectation 
values of the traced spatial and time-like link matrices in the Landau gauge ||. The 
actual values are established self-consistently by an appropriate iteration procedure. 
The choice of Landau mean-link tadpoles is motivated by evidence that their use leads 
to better continuum behaviour than can be obtained from plaquette tadpoles [||, 0, g|. 

The renormalized anisotropy is \R = a s/&t where a s and a% are the spatial and 
time-like lattice spacings respectively. The value of \R can De obtained in several 
different ways. In a separate paper we have determined it by fitting torelon and glueball 
dispersion relations measured on the lattice and by measuring the heavy quark potential 
in both the fine and coarse directions ||. In this paper we use the bare parameters = 
1.8 and x = 6 together with the appropriately tuned tadpole parameters u s = 0.7216 
and ut = 0.99208, with corresponding value \R = 5.32(2) . From the T(1P— IS") mass 
gap measured on the lattice we find a t _1 = 4503(46) MeV and a' 1 = 819(9) MeV. 

Our calculations are a continuation of previous work || |l0[ in which we investigated 
the effect of including all terms in the quark NRQCD Hamiltonian up to 0(mw 6 ) . 
However in measuring the heavy hybrids we have found that terms of 0(mv 4 ) and 
higher have negligible effect on our results. We therefore drop such terms in order to 
reduce the computational load. This circumstance is consistent with the observation 
that the heavy quarks move rather slowly in the hybrid state. 

The contribution from the first spin dependent term is large however and central 
to our results. Accordingly we retain this term and the 0(a s ) improvement to D 2 , and 
use the quark Hamiltonian H = Hq + 5H, where 

a, = ° 2 



A 4 g 

, Q _ C\ 

8m^ 2ml) 



5H = -c --ci^^a.B + Oimv*). (1.2) 
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The -B-field is constructed from an improved form for F^ u [11|. We take the tree-level 
value Co = 1 . 

The thrust of our investigation is to study the dependence of the hybrid mass 
spectrum on c\. This parameter may be subject to a large renormalization that forces 
its departure from the tree level value c\ = 1 . The states considered are the low-lying 
magnetic hybrids which are classified as the 1 , which has quark spin s = 0, and 
the , 1 , 2 ^ , which have s = 1. The 1 is non-exotic and is the only one of 
these states that can be produced directly in e + e~ collisions. Of the others, the 1 h is 
exotic. 

In section [2] we outline the NRQCD approach to calculating heavy hybrid states, 
and discuss the operators and smearing used. In section [5| we discuss finite-size effects 
and present the results from our simulation. In section ^ we discuss the bag model cal- 
culation of spin effects and inter-state mixing in the Born-Oppenheimer approximation. 
In section |7| we present our discussion and draw our conclusions. 

2 The NRQCD approach 

The NRQCD approach to computing the heavy hybrid propagator is well known The 
Euclidean time evolution equation for the quark propagator G(x, t; y) is computed 
according to the prescription: 

G(x, t + 1; y) = (l - °^f\ n C4(x) (l - " (1 - a t 5H)G(x, t;y) t>0, 



where Hq and 5H are given in eq (|1.2|) and n = 2 for our simulations. The quark 
propagator is then determined by its initial distribution, G(x, t = 0; y) = T(x, y) . The 
choice of initial distribution is governed by the use to be made of the propagator. We 
use propagators computed from the following in initial distributions: 

r L (x, y ) = 
r 5 (x, y ) = 

Tf(x,y) = 

Tf(x,y) = Z^l + ^-j ^ 3 )(x,y), (2.1) 

where Di is a (discrete) covariant derivative. The propagator resulting from the local 
source, r , is combined with the each of the other smeared propagators, after hermitian 
conjugation, and traced with appropriate spin matrices to create the correlators for the 
set of operators listed in Table [j]. 

The correlators appropriate to hybrid states are computed from Yf . An analysis 
of hybrid states in the Born-Oppenheimer approximation indicates that the quark and 
anti-quark tend to remain apart from one another because they experience an angular 
momentum barrier and a color repulsion (the qq system is in a color octet). The "donut" 
smearing incorporated in F H through the application of a final D 2 , was successfully 
invoked in order to improve the overlap of the operators with the widely separated qq 
system in the hybrid. 
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The forms of the operators we use to interpolate the magnetic hybrid states are 
shown in Table |j. 
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Table 1: The forms of the operators used to interpolate the magnetic hybrid states we study 
and their J classification. Here I and s are the QQ angular momentum and spin, and j is 
the gluon angular momentum. 



The charge conjugation (C) and parity (P) of the quarkonia and hybrids, satisfy 

pure magnetic electric 

quarkonium hybrids hybrids 

p= (-i) m (-i)^' +i 

where I and s are the qq angular momentum and spin and j is the gluon angular 
momentum. The S,P and some of the hybrids are listed in Table 0. 



3 Computational details 

In the source and sink operators (Table |l]) all links were smeared using covariant 
smoothing in a manner similar to that described in eq( [2.lD with an appropriately de- 
fined D 2 . The B fields were computed from the smeared links and were themselves 
subsequently smeared in a similar manner. The range for the smearings was typically 
I = 2 on gauge links and I = 3 on the quark source. The effectiveness of the smearing 
was confirmed by a reduction in statistical errors by a factor of three or more. To reduce 
statistical errors to within 1% we generally required between 40,000 and 60,000 inde- 
pendent correlators. We used a 50 sweep Cabbibo-Marinari thermalisation to generate 
the initial gauge fields, and then employed 10 sweeps between measurements. 

The computations were performed on the Hitachi SR2201 of the Cambridge High 
Performance Computing Facility and at University of Tokyo Computer Centre. They 
were farmed across 64 processors each of which returned statistically independent re- 
sults that were used to construct the full covariance matrix. In the standard way we 
performed a correlated fit to the correlation functions of the relevant operators in order 
to determine the masses of the quarkonium and hybrid excitations. We confirmed that 
10 sweeps between measurements was sufficient for decorrelation. 
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The timings for the calculation was such that for the hybrids an 8 hour job produces 
about 80 correlators per processor for the 6 3 x 36 lattice, about 25 for the 8 3 x 48 and 
5 for the 10 3 x 60. 



4 Finite size effects 

We first examined the sensitivity of the calculation to finite size effects. In Table || we 
show the mass gap from the T(15) to the 1 hybrid (without spin effects) measured 
by us on lattices L 3 xT for a range of values of L . Only when L is reduced down to 4 
is there any detectable effect from the finite spatial lattice volume. This is consistent 
with the results of Manke et al || . 

since a hybrid is typically 2-3 GeV" 1 g| across, its size is ~ Aa s . For L = 6 it 



covers a substantial part of the lattice volume. It is surprising therefore that finite size 
effects are so small since we might expect that a quark would react more strongly to 
the nearby image of an anti-quark across the boundary of the spatial lattice than to the 
more distant anti-quark within the the spatial volume. However a strong effect depends 
on a transition in which the gauge field flux tube switches from the original quark 
anti-quark linkage to one crossing the spatial boundary. The corresponding transition 
amplitude involves the expectation value of a closed Polyakov loop that encircles the 
torus. This vanishes in the confining phase because of the Z3 centre symmetry of the 
action. We argue therefore that our results are consistent with a picture in which the 
major interaction across the boundary of the spatial box is a weak Van der Waals-type 
force. 

There are other re-linking transitions, in the case of SU(3), involving the multiple 
wrapping of flux lines round the torus, that do not vanish. However we expect them 
to be suppressed because of the high energy of the associated flux tubes. 

For the rest of our investigation we worked with L = 6 as the most economical 
choice of lattice consistent with negligible finite size effects. 



5 Results 

We use as an energy baseline the spin averaged T(15) state. Table || shows our results 
for the dependence of this baseline on c\ . They are consistent, with quadratic behaviour 
in ci reflecting the fact that the shift is due to a quark mass renormalization. 

Since the (spin averaged) IS — IP mass difference AMps in lattice units is known 
to be insensitive to the quark mass mi, we use it to determine the temporal lattice 
spacing in physical units. From the assignment AAlps = 440 MeV we infer that 
a^ 1 = 4503(46) MeV . From other work [|| we know that xr = 5.32(2) yielding 
a^ 1 = 819(9) MeV . We determine nib by measuring the mass of the T(1 3 S') state from 
its dispersion relation and then requiring it to be consistent with the mass measured 
directly. We find m b = 4.75 GeV . 

Table || shows our results for the masses of the hybrid states. We also show here 
the behaviour of the mass with variation of the c\ co-efficient in 5H, which controls the 
magnitude of the <r.B term. These results are illustrated in Fig. 

For each choice of c\ (defined in eq (|l.2|) ) we compute the s = 0, 1, / = quarko- 
nium correlators and the hybrid correlators for the states listed in Table Although, 
in principle, we could consider a matrix of source-sink operator pairs and fit the cor- 
responding matrix correlator, we found it sufficient to consider the one correlator for 
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(L,T) 
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(6,36) 


(8,48) 


(10,60) 


{M H - M r )a t 


0.281(7) 


0.244(7) 


0.244(9) 


0.247(9) 



Table 2: Dependence of the mass, Mjj, of the 1 hybrid state on the finite size, L. In all 
cases ci = 0, a^ 1 = 4503(46) MeV. 
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Table 3: Mass of the 15*3, IS® and the spin-averaged value for different values of c\. 

which both source and sink were smoothed as described above. We found that the 
S-wave quarkonium correlators were well fitted with two-exponentials over the range 
(in lattice units) 10 < t < 50. The hybrid states are much heavier and the signal for 
the propagator was too noisy for t > 15 but we obtained excellent fits with a single 
exponential for typically t > 5. An representative effective mass plot for the 1 — h for 
ci = — 1 is shown in Fig. ^ which has a plateau extending from t = 4 — 10. The dashed 
lines shown the upper and lower error bounds from the fit. 

The 1 hybrid is of particular interest since it will appear as an intermediate 
state in e~e + scattering at .B-meson factories. It is noteworthy that its excitation 
energy above the ground state is unchanged by the spin interaction independently of 
its strength. The absence of a term 0{c\) implies that there is no direct energy shift 
induced by the spin term in the Hamiltonian and the absence of a term 0(cf ) implies 
that there is no mixing with quarkonia even though this is possible on the basis of 
quantum numbers. We explain below why we feel that this result is compatible with 
the bag model of hybrid states. Our prediction for its mass is 1644(17) MeV above the 
ground state. This is compatible with, if somewhat higher than, the value from some 
other studies || 10, 13, pf| . It would be of great interest to observe experimentally a 



state at this energy level. 

In contrast the masses of the spin one states are strongly affected by the spin term in 
the NRQCD Hamiltonian. Of these the 1 ^ is exotic and cannot mix with quarkonium 
states. For this state we have investigated the dependence of its mass on c\ for both 
positive (physical) and negative (unphysical) values. This allows us to disentangle 
contributions that are linear and quadratic in c\. From Table || we see that the effect 
of the spin term is large. For the tree level value c\ = 1 the 1 h is depressed by ~ 130 
MeV. A large renormalization of c\ because of radiative corrections means that ci = 2 

is a possibility, a value which causes the 1 |_ to be depressed by almost 300 MeV to 

1265(26) MeV above the spin-averaged 15 QQ state. A simple fit to the 1 h masses is 

a t {E H - E q q) = 0.369(4) - 0.026(2)ci - 0.09(l)c? , (5.1) 

with x 2 = 0.3. These effects are unexpectedly large since spin effects depend on the 
quark velocity which is low in hybrid states. 

The non-exotic ^ and 2 h states can mix with quarkonium states and, in prin- 
ciple, there can be contamination of the signal which would make these states hard to 
pick out in the simulation. However, as in the case of the 1 we find no such problem. 
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Table 4: Mass of the hybrid states with variation of the spin coupling magnitude, accurate to 
0{mv 6 ) with [3 = 1.8, am b = 4.75, aj 1 = 4.503(46) GeV. 

On the basis of a bag model calculation this mixing is indeed found to be very small 
and certainly at a level of less than 1% in our simulations. The h is depressed even 
more than the 1 ^ and for c\ = 2 is only ~ 1000 MeV above the QQ ground state. 
The 2 ^ rises for c\ > although the magnitude of the shift is not as great as that 
for the other states. This is compatible with a change of sign of the coefficient of c\ in 
the fit eq( |5,l| ) which causes this term to compete with the quadratic contribution. We 
shall argue in the next section that at least a part of the coefficient of the term linear 
in ci is proportional to L ■ s which does indeed have this qualitative behaviour. We 
have not investigated the 2 for c\ < 0. 



6 Born-Oppenheimer Model 



Hasenfratz et. al [|12j l and Kuti and Morningstar [14] have stressed the usefulness of a 
Born-Oppenheimer picture of hybrid states in the bag model of hadrons. We suggest 
that our results can be consistently interpreted in this approach. We investigate the 
effect on mass shifts of a gluon exchanged between the quark and the anti-quark, and 
the mixing effect induced by transitions to a two-gluon state. We find the former to be 
very small while the latter can provide a sizeable contribution that goes some way to 
explaining our observations on the mass shifts of the hybrid states. We also estimate 
the mixing of the hybrid and quarkonia states and find this to be negligibly small for 
all reasonable values of c\. 

For simplicity we assume that the bag is spherical. The wavefunction for a hybrid 
of total orbital angular momentum (L,m) in the BO picture is 

*Lm(ry 9 ) = f^(r)D^ A (n)Aj A (r' g ) , (6.1) 

where r is the quark-antiquark separation vector with orientation f2 = (0, (f>), r' g is the 
gluon position vector in the body frame, j is the angular momentum of the gluon mode 
with projection A onto the body-fixed z-axis, and A = 0, 1 label TM and TE modes, 
respectively. The total angular momentum is J = L + s. States of definite parity are 
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given by 

= ^Ki + rfj • (6-2) 
The parity (P) and charge-conjugation (C) are given by 

p = n(-l) L+x+1 , C = r)(-l) L+j+s+1 . (6.3) 

Note, that these assignments do not agree with those quoted in the original reference 



[12], but they do agree with Le Yaouanc et. al [15] and Rakow [|16| who remarked on the 
discrepancy. We have carefully checked that these formulae are correct. The important 
point is that the PC classification is independent of A. States with |A| = 0,1,2,... 
are labelled respectively by E, II, A, ... . The lowest lying states include those with 
L = l,j = 1 on which we shall concentrate. There are both E and II states with the 
same J PC quantum numbers, namely, 1 with s = and h , 1 ^, 2 ^ with s = 1. 

There are two simple processes resulting from the a.B interaction, gluon exchange 
between the quarks, and gluon exchange between the hybrid gluon and one or other of 
the quarks. The former is a weak effect because it involves the behaviour of the QQ 
wavefunction near the origin where it is suppressed both by the angular momentum 
barrier and a repulsive potential due to the octet nature of the QQ state. An estimate 
yields a result of order 1 MeV or less. 

The large effect is due to gluon exchange between the constituent gluon and the 
quarks. This effect is large for two reasons: the colour factors are large and it depends 
on the quark wavefunction in the bulk and not on its value at the origin. We outline the 
calculation below and follow the formalism of Close and Horgan [jlT], [0| and Barnes, 
Close and Monaghan [jE], [2^] . The latter authors present a calculation of this spin effect 
in the MIT bag model and we do not reproduce all the details here. The calculation 
is for the contribution of second-order perturbation theory to the mass shift. The 
intermediate states are summed over and, in our calculation, are BO states with two 
constituent gluons. We examine the contribution from one intermediate state with 
jl = j2 = 1, Ai = A2 = 1. the diagram for this process is shown in Fig. |3| and 
corresponds to the exchange of a j = 1 electric gluon. In this case, it turns out that 
L' = L, s' = s and Ai + A2 = A, where A labels the initial state gluon and Ai, A2 label 
the intermediate state gluons as shown in Fig. ||. Since AMh also does not depend on 
the sign of A we consider A = 0, 1 which correspond to the E and II states, respectively. 
The bag is taken to be spherical of radius Rb- The interaction Hamiltonian and the 



expressions for the gluon modes, A^, and associated B fields are given in [19]. The 



field operators contributing to the mass shift are expanded on this mode basis using 
the usual creation and annihilation operators. The two-gluon intermediate state can 
be written in the BO approximation in a form similar to that in eqn. (|6.1|). After some 
manipulations of the D-matrices involved and some algebra the mass shift takes the 
form: 

AMh = (-1) A e A ■ (e* Al Aeyi.s, (6.4) 

m b K B 

where the eA=o = (0,0,1) and e A= ±\ = (— 1, =pz, 0)/\/2 and C is an overall constant 
which includes colour factors, the characteristic bag constants and the radial overlap 
integrals for the quark wavefunctions. The bag radius is Rb ■ 

Our estimate of C along the lines of |20[], yields a relatively large value C ~ 1.13 . 



The factor e\ ■ (e\ Ae! ) is typical of the three-gluon vertex and is discussed in IS, 20] 
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One important result is that AMh = for A = 0, i.e., for the E state. The outcome 
can be summarized as: 



AM H = Kcia s x 
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Our calculation is, by its very nature, approximate. For example, the bag is not spher- 
ical but elliptical, a correction included by Kuti and Morningstar jnj] when calculating 
the QQ BO potentials. Including this correction will increase the effect. Also, the bag 
radius is not fixed but varies with the quark motion in the BO approximation. The 
value we adopt is therefore associated with the most likely quark separation and we take 
R B ~ 1.5 - 2 GeV -1 . With m b ~ 4.75 GeV we find K ~ 60 - 110 MeV. With c x = 1 
and a s = 0.4 we find the scale of the induced spin-orbit splitting is ~ 25 — 40 MeV to 
be compared with ~ 117 MeV for the linear term in c\ from eq(|5.1|). The value chosen 
for a s is typical of bag calculations and is consistent with av(/J>) at a suitably chosen 
scale based on the B-quark mass |2l| . While our prediction is too small our calculation 
is clearly incomplete. There are many other intermediate states that can contribute 
with similar magnitude of shift. The contributions of diagrams involving the interac- 
tion of the quarks with the constituent gluon are large partly because the colour factors 
are large and partly because there are many such terms. At 0(as) there is a similar 
diagram to the one analyzed with j = 2 electric gluon exchange. At 0{a 2 s ) there are 
diagrams involving the gluon four-point coupling, such as shown in Fig. which has 
many orderings and gives contributions linear in c\. There are many similar and other 
higher-order diagrams which contribute to the c\ term. It cannot be excluded that 
there are also non-perturbative effects as well. The dependence of the induced mass 
shifts on the quantum numbers will be generally more complex than the L ■ s derived 
explicitly above, which nevertheless gives a qualitative picture of the pattern observed. 
We are therefore confident that our calculation sets a lower bound on the effect. 

In the previous section we remarked that the 2 *~ rises above the c\ = reference 
mass by an amount smaller in magnitude than the shift of either the h or 1 h below 
this point. This could be due to a cancellation between the linear and quadratic terms 
in eq( |5.1| ). However, our operators do not distinguish between the E and II states since 
they project out definite quark angular momentum Iqq = which corresponds to an 



equally weighted mixture of A = 1,0,-1. Based on the potentials of Kuti et. al [14] 
we see that the E (which has the same PC assignments as the n) is about 150 MeV 
above the II. Since the E is unaffected by the spin-orbit term it may be that we are 
seeing the 2 ^ E state and not the II state because this latter has risen to a higher 
mass. We are testing this suggestion by constructing operators diagonal in A and the 
results will be presented in a subsequent publication. 

In order to confirm our observation that mixing with quarkonium states is negligible 
we have estimated the mixing in the BO approximation in the bag model. Since both 
source and sink contain a gluon field the gluon must be annihilated and subsequently 
recreated. A diagram of the kind that is relevant is shown in Fig. ||. The states that 
mix must have opposite spin and so the mixing is by a spin-flip operator, in this case 
it is proportional to (<tq — ctq). We calculate this term for the s = 1, ground-state, 



9 



hybrids mixing with the s = quarkonium state. A term of this kind has been discussed 
earlier by Ono [22] and Le Yaouanc et. al [15]. Our result is compatible with theirs 
and we find that the mixing amplitude from this diagram is ~ 10" 4 . Ono [E2| gives 
~ 2 x 10 -4 . The hybrid-QQ mass difference in lattice units is ~ 0.36 and so the effect 
of mixing at t = 10 is ~ 10 -4 x e~ 3 ' 6 ~ .4%. It is unlikely that mixing because of 
excited hybrids of the same quantum numbers will contribute appreciably and so we 
expect that mixing has negligible effect on the signal. 

We believe that we have demonstrated that spin effects in heavy hybrid mesons are 
large and that the explanation is that the states expanded in the Fock basis contain 
more than one gluon constituent with an appreciable probability. Our partial estimate 
based on one intermediate state is ~ 20 — 30% of the observed effect and shows that 
the order of magnitude can be explained. We are in no doubt that the complete sum 
over intermediate states will contribute to the shortfall and, because the effect at lowest 
non-zero order is so large, it is sensible to suppose that higher order corrections are 
appreciable. The most easily observed 1 state is unaffected since it has s = but 
indirect methods of production should eventually detect the s = 1 states and then our 
observations can be tested. 



7 Discussion and Conclusions 

Ours is the first study on an anisotropic lattice that incorporates tuned Landau tadpoles 
in the action. The choice of Landau tadpoles was motivated by the observation that 
they yield better scaling behaviour for observables than plaquette tadpoles ||. 

We have confirmed that the hybrids are at a levels in the range 900 - 1750 MeV 
above the IS spin averaged ground state. This is consistent with our previous work 
and, allowing for systematic uncertainties in at inevitable in a quenched calculation, is 
consistent with other evaluations. 

The important new result of our calculation is that the effect of the spin interaction 
on the hybrid spectrum is unexpectedly large. We measure energy level splittings 
~ 100 MeV and greater. The results are exhibited in Table |] for a range of values of c\ 
including, in some cases, unphysical negative values that help to establish parameters 
for the analytic form of the dependence on this coefficient. 

In the absence of the spin interaction (ci = 0), all the states we measure are 
degenerate. Because the exotic cannot mix with quarkonia, neither do the other states 
in the multiplet. When the spin interaction is switched on two effects can occur, 
(i) a shift of the energy levels of the hybrids and (ii) a mixing of the non-exotics 
with appropriate quarkonium states. We claim that our results are consistent with a 
scenario in which there is a direct effect on most of the energy levels but little mixing 
with quarkonia. 

At tree level with no lattice artifact corrections, c\ = 1. In deriving the NRQCD 
approximation from an interacting gauge theory we would expect corrections to c\ of 
0(a s (ma s ) 2 ). These could be large enough to yield c% ~ 2 . Hence, the spin effect 
we have simulated here may be a reasonable prediction even for c\ ~ 2 . Because 
the terms in the NRQCD Hamiltonian which are 0(mv 4 ) and higher have neligible 
effect on the hybrid levels we can expect the effect of radiative corrections to be limited 
to a renormalization of c\ which gives hope that the spin effects observed here are 
quantitatively correct. 

It would be of great interest if any of these hybrid states could be observed although, 
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except in the case of the 1 which can be directly produced in e + e~ collisions, this 
would be experimentally difficult. 

Future work will concentrate on using source and sink operators which have the 
quantum numbers of the Born-Oppenheimer approximation. This will allow us in 
particular to disentangle the II and £ states. 
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Figure 1: The variation of the hybrid masses with C\ for the 1 ,0 + , 1 + , 2 
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Figure 2: Effective mass plot for the 1 + hybrid with ci = — 1. The dashed lines show the 
upper and lower bounds of the fit. The plateau extends from t=4-10. 
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Figure 3: A representative graph for the second-order perturbation theory calculation discussed 
in section ^ for the mass shift due to the inclusion of the cr.B term in the NRQCD Hamiltonian 
(1.2). This spin operator is represented by the filled box and terms in the spin independent 



part of the interaction Ham iltonian are represented by filled circles. This contribution is to the 
linear term in c\ in eq(fTl|) and is 0(0:5). 




Figure 4: A higher-order, 0(a|), perturbation contribution to the linear term in ci in eq( |5.l| ). 
Symbolism is the same as in Fig. 0. 




Figure 5: A representative graph for the mixing contribution to the hybrid- hybrid correlator 
from quarkonium intermediate states. 
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